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Abstract 

The solutions of the Heavenly Equations and their generalizations are considered. 

1 Method of solution 

To integrate the partial nonlinear first order differential equation 



' Hi ^1 "^1 fxi fyi fzi fwi fxxi fxyi fxzi fyyi ••) (1) 



can be applied a following approach [1]. 

We use the parametric presentation of the functions and variables 



f{x,y,z,w)-*u{x,t,z,w), y ^ v{x,t,z,w), fx^Ur,-—v^, 
X ; ^* 

fw ^ ^10) fz ^ "^2 l^z^ fy ; (2) 

Vt Vt Vt 

where variable t is considered as parameter. 
Remark that conditions of the type 

fxt ftx J fxz fzx 5 fxw fwx • • • 

are fulfilled at such type of presentation. 

In result instead of equation (1) one get the relation between the new variables u{x,t, z,w) 
and v{x, t, z, w) and their partial derivatives 

V, Ux, Uz, U^, Ut, Vx, Vz, Vnj, Vt...) = 0. (3) 

In some cases the solution of such type of indefinite equation is more simple problem than the 
solution of the equation (1). 

We apply this method to integrating of the equations having applications in various branches 
of modern mathematical physics. 
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2 The first heavenly equation 

The Plebanski first heavenly equation 

arise in connection with the 4-dimensional spaces endowed with the metrics 

ds^ = 29xydxdy + 29zwdzdw + 29xwdxdw + 29zydzdy (5) 

which are simultaneously Ricci flat and self-dual ([2]) 
In result of applying of the transformations 



and 



9{x,y, z,w) = u{x,t, z,w), y = v{x,t, z,w) (6) 

d9{x,y,z,w) d (f^u{x,t,z,w))£v{x,t,z,w) 
ax ax ^v{x,t,z,w) 

de{.r. y. z,w) d . , (|»(:r. f . - . »0) -|r(:r. f . - . »■) 

= —u{x,t,z,w) - ^ Q ' ^ , 

az oz ^v{x,t, z,w) 

d9{x,y,z,w) d . , (f^u{x,t,z,w))£v{x,t,z,w) 

a ^ ^^(^' ^' ^' - ^ a , \ ^ ' 

aw aw ■^v{x,t, z,w) 



d9{x,y,z,w) ^u{x,t,z, 



w 



dt ■§-^v{x,t,z,w) 

with corresponding expressions for the second order derivatives at the equation (4) one gets the 
relation between the functions u{x, t, z, w), v{x, t, z, w) and their partial derivatives 

F{u,v,u^,ut,u^,u^,...) ^0. (8) 

The corresponding explicit expression for the function F can be obtained with the help of the 
MAPLE and we omit its in view of its inconvenience. 

To cite as an example some of the simplest reduction of full expression (8). 
The substitution ^ 

v{x, t, z, w) — t—u!{x, t, z, w) — uj{x, t, z, w), 



u[x^ t, Z, W) = —UJ[X, t, z, w) 



give us the equation 



/ 92 \ ^ Q2 
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— z, w) 1 -Q^'^ix, t, z, w) = 0. (9) 



It has the particular solution 

u;{x, t,z,w) = z + e^+'^t (_Ci sinh(t-^) + _C2 cosh(t-i 

where _C1 and _C2 are the parameters. 

After substitution of this expression into the relations 

d d 
y — t—u!{x, t, z, w) + u!{x, t, z, w) = 0, 9{x, y, z, w) — —u!{x, t, z,w) = 

we find the parametric presentation of variables (in particular case _C2 = 0) 

y + e^+^'.Ci sinh(t-^) + z = 

and 

0{x, y, z, w)t - e^+^_Ci t cosh(r ^) + e^+^_Ci sinh(r ^) = 0. 

Elimination of the parameter t from both relations give us the explicit expression for the 
function 



p-2x-2wii2 I 9 p-2x-2w,,iy _i_ p-2x-2w -,2 _i_ /7/2 

e{x, y, w) = _C1 ^ 

[y + z) arcsmhy— 



_C1 

which is solution of the first heavenly equation (4) . 
The substitution 

d 

u{x, t, z, w) — t—uj{x, t, z, w) — uj{x, t, z, w), 

I \ 9 . , 

V{X, t, Z, W) — —U}[X, t, Z, W) 

into the expression (8) lead to the equation 

-Q^^i^A^^^'^) (10) 

It has solution in form 

u!{x, t, z, w) = A{x, t, z)w 

where 

-—A{x,t,z)^Q. 
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Prom here we find a simplest solution 

A{x,t,z) = ln(- -l) + x. 
z 

In result we get the relations 

yz + w 

i — 



and 

2; f z \ 

—t + ln( )t — ln( )z + xt — xzj {t — z)~ . 

Elimination of parameter t from these relations give us the explicit expression for the function 

w 

9{x, y, z,w) — yz + w — w ln( — ) — wx 

yz 

which is solution of the first heavenly equation (4) . 

By analogy can be constructed a more comphcated solutions of the equation (4) . 

3 The second heavenly equation 

The Plebanski second heavenly equation has the form 

-0{x, y, z, w) + ^-^0{x, y, z,w)+ [ ^0{x, y, z, w) Tr^O{x, y, z, w)- 



awax ayaz \ox^ ) oy^ 

= (11) 

It describes the properties of the Ricci flat a 4-dimensional spaces endowed with the metrics 

ds^ = —Oxxdz^ + 26xydwdz — OyydnP' + dwdx + dzdy (12) 

satisfying the conditions of self-duality for corresponding Riemann tensor Rijki of the space [3]. 
To integrate the equation (11) we transform its in another form with accordance of the rules 

(2). 

In result of such transformation one gets the relation between the functions u{x, t, z, w) and 
v{x, t, z, w) and their partial derivatives 

F{u, V, Ux, Ut, U^, Uyj, Uxx, Uxt, Uxz--) = 0. (13) 

The next step is receiving the p.d.e from this relation using the substitutions similar to 

d d 

u{x, t, Z,w) — t —U!{x, t, Z, w) — U!{x, t, Z, w), v{x, t, Z, w) — —U!{x, t, Z, w), 

or 

d d 

v{x, t, z,w) —t —uj{x, t, z, w) — uj{x, t, z, w), u{x, t, z, w) — —uj{x, t, z, w), 

or the more complicated. 



4 



As example we get from the relation (13) the equation 

/ Q2 \ Q2 Q2 Q2 

with the help of the first substitution. 

The simplest solution of this equation is 

(jJ\X^ t, z, w) — h z 

X 

and corresponds the solution of the equation (11) 

9{x, y, z, w) = -1/4 y ^^"^^^ 

w 

after the return at the assumed function 0{x, y, z, w). 
The substitution 

U!{x, t, z, w) — B{t, z, w) + xC{t, z, w) 
into the equation (14) give us the system 

Its simplest solution is 

C[t,z,w) — — 



and 

B 



{t, z,w)^w (^_F2{z) + j _F1 {z (l + e~*))eHt + j _F1 {z (l + e'*))dt^ , 



where _Fi{z) are arbitrary. 

Using these expressions we get the function u!{x, t, z, w) from (14). 
In particular case we have 



/ / N / /\ x (e* + 1) w 

uj{x,t,z,w) = {_F2{z) +ze^ + 2zt- ze~^)w + — ^ —. 

^ ^ z 

Using this expression and the relations 

d 

y - -Q^^i^^t^z^w) = 0, 
d 

9(x, y, z, w) — t—uj(x, t, z, w) + u(x, t, z,w) — 
at 

we find 

yz — wz^e* — 2 wz^ — wz^e~* — wxe^ — 

and 

/ f f f A V V wx — twxe* + wxe^ 

[-twe-twe * + we-we ^) z + 9{x,y, z,w) + W-F2{z) -\ = 0. 

^ ' z 

Elimination of the parameter t from these relations give us the function 9{x,y, z,w) which is 

solution of the second heavenly equation (11). 
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4 The Dunajski heavenly equation 

The Dunajski heavenly equation has the form of the system 

Q2 Q2 

-f{x, y, z, w) + ——f{x, y, z, w) + 



awox oyaz 



-L(a;, y, w) + ^-^L(x, ^, w) + ^f{x, y, z, w) ^L{x, y, z, w)+ 



awox oyaz \ox'^ j ay 

+ ^)) ^^(^' ^) - 2 [-id^fi^^ ^)) ^) = 0- (15) 

It describes the properties of a four-dimensional Riemannian spaces having the metrics 

ds"^ — —fxxdz^ + ^fxydwdz — fyyduP' + dwdx + dzdy (16) 
with condition on the Ricci-tensor of the metric (16) 

where Tij is a Maxwell stress-energy tensor of electromagnetic field [4]. 

The conditions of self-duality for the Riemann tensor Rijki of these spaces are equivalent the 
equation (15). 

In connection with this it is interested to note that the equation (15) coincides with condition 

(^1212 = 0, 

for the component C1212 of the Weyl tensor of the metrics (16). 

To construction of the solutions of the equation (15) we transform its in accordance with the 
rules (2). 

In result we get the relation 

^{U, V, Ux, Ut, U^, Uyj, Uxx: Uxu Uxz---) = (17) 

which is reduced to the p.d.e.'s when between the u and v exist functional dependence. 
Let us consider an examples. 
In the case 

d d 
u{x, t, z, w) = t—u!{x, t, z, w) — u{x, t, z, w), v{x, t, z, w) = —uj{x, t, z, w) 

one get the p.d.e. on the function uj{x, t, z, w) 

M{UJ, UJ^, UJt, UJz, ^w, •••) = 0. 

Its simplest solution can be presented in form 

ijj{x^ t, z, w) — A{t) + xzt + w, 
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where the function A{t) is defined by the equation 



In result we get the solution 

u;{x, t, z, w) = 2/3 — ^ '- V-C3t + + xzt 

_0-Z 

from which is followed 

t = 

-_C2 - 1/2 _C1 + _Ci y_C3 + _C1 yxz - 1/2 _C1 _C3^ - _C1 _C3 xz - \/2_Cl xV, 

and 

f{x,y,z,w) = 



-t^V2 + tV2-C2 + 2 V2-C2^ + 3 ^f^IcTitV^C^ + 3 w^-_Cl {t + _C2) 

Elimination of the parameter t from these relations give us the function 

f{x,y,z,w) = l/Q _C1 x^z^ + (-1/2 _C1 x^y + 1/2 _C3 _C1 x^) z^+ 

+ ((l/2 _C1 _C5^ + _C^) X + 1/2 x_Cl y'^ - _C1 _C3 xy) z - 1/6 y^.Ci + 1/2 _C3 _C1 

-w + (-1/2 _C1 _C3^ - _C2) y + 1/6 _Ci _C3^ + _C3 _C2 - _C4 

satisfying the equation (15). 
The substitution 

d d 

v{x, t, Z, w) = t—U!{x, t, Z, w) — Lj{x, t, Z, w), u{x, t, Z, w) = —Uj{x, t, Z, w) 

lead to the p.d.e. on the function uj{x, t, z, w) 

N{uj,UJ^,UJt,UJz,l^w, •••) = 0. 

having the solution 

u!{x, t, z, w) — Ae^w -\- {x -\- z) t. 
This solution give us the following presentation 

y - tAe^w + Ae^w = 

and 

f{x, y, z, w) — Ae^w — x — z = 
Elimination of the parameter t from these relations give us the function 

/(x, 2/, z, w) = e"'='^^(^)+^/lw + X + Z 
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or 

f(x,y,z,w) = 

y + xLambertW{— — ) + zLambertW{— — )j {LambertW{— — )j , 

satisfying the equation (15). 

The LambertW function in these formulas is defined by the equation 

LambertW {x)exp{LambertW (x)) ~ x. 

From here it is apparent how to construct another examples of solutions of the equation (15). 

5 Six-dimensional generalization of the first heavenly equa- 
tion 

The first heavenly equation in dimension Dim=6 has the form 

— -— Afx) A(x) A(x) — A(x) A(x) A(x) I 

\dudz J \dvdw J dxdy \dudz J \dwdx J dvdy 



^(^) ^-^^(^) ^-^^(^) - ^-^^(^) ^-^^(^) ^-^^(^)- 



^dydz J \dwdx J dudv \dwdz J \dudv J dxdy 



-A^) ^Mx) ^A{x) + ^Mx) ^A{x) ^Mx)- 



\dydz J \dudx J dvdw \dudx J \dvdy J dwdz 

-1 = 0, (18) 

where A{x) — A{x, y, z, u, v, w). 

It arises as condition on the six-dimensional metrics 

ds^ = 2 (t^^A(x)] dxdy + 2 ( tt^Mx)] dxdu + 2 ( ^ Mx)] dx dw+ 
\oxdy J \ouox J ydwox J 



to be a Ricci-fiat 



Rij — 0. 



The transformation of the equation (18) at the new form in accordance with the rules (2) lead 
to the relation 

^{U, V, U,, Ut, U„ U^, U^, U,t, U,z..., V,, Vt, V„ V^, V^, V,,, V,t, y.z-) = 0. (19) 
This relation moves to the p.d.e. 

/ 92 \ f \ 92 / 92 \ / 52 \ 92 

- is^'^'^V l^^'^'V ^"^"^ " v^"^"V 9^"^"^^ 



+ i^^'^^^V la^'^'^j a^"^"^ " ia^'^'^V l^"^"^J a^"^"^" 
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at the substitutions 



d 

U (x, t, z, u, V, w) = t—Lulx, t, z, u, V, w) — Lulx, t, Z, U, V, w), 
at 

d 

V{x, t, z, u, V, w) = -^i^{x, t, Z, U, V, w), 



where oo{t) = uj{x, t, z, u, v, w). 

The solutions of last equation allow us to construct the solutions of initial equation (18) 

Let us consider an examples. 

The solution of the equation (20) of the form 

uj{x, t, z, u, V, w) — Bxt + uv + zw + 1/2 Bt^ 

corresponds the solution of equation (18) 

, , N . -y^ + 2 yBx - B^x^ + 2 uvB + 2 zwB 
A{x, y, z, t, u, V, w) = -1/2 . 

The solution of the form 

/ , \ 1 {xz + 1) 

uj[x^ t, z, u, V, w) = ln[ h _L2)u + vw 

z 

lead to the solution of equation (18) 

Aix, y, z, t, u, v,w)^ - ly-Cl xz + y^CS z- -Clu + ln( "^^ ^ )u-Cl + vw^Cl \ ^Cl 

V yz J 

The solution of the equation (20) of the form 

co{x, t, z, u, V, w) = sinh(a; + t) {zu + vw) 
give us the solution of equation (18) 



. , , / —y + ZU + vw y + ZU + vw 

A{x, y, z, t, u, V, w) = -W ■ W ■ vw - yx+ 

V ZU + vw V zu + vw 



, , y , / —y + zu + vw y + zu + vw 

+yarccosn{ ) — W W zu. 

zu + vw V ZU + vw V ZU + vw 
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6 Six-dimensional generalization of the second heavenly 
equation 

The second heavenly equation in dimension Dim=6 can be obtained as generahzation of the metric 
(12). 

It looks as 

ds'^ = A{x)dx'^ + 2 B{x)dxdy + C{x)dy^ + 2 E{x)dxdz+ 

+F{x)dz'^ + 2 H(x)dzdy + dxdu + dydv + dzdw. (21) 

The Ricci tensor of such metrics has a fifteen components among which are the nine components 
expressed trough the values 

where x — {x, y, z,u,v,w). 

Equating these expressions to zero we get the equations for determination of the components 
of the metrics (21). 

d d d 

—£"(0;, y, z, u, V, w) + g;^H{x, y, z, u, v, w) + V, ^, v, w) = 0, 

d d d 

— 5(x, y, z, u, V, w) + g:^C{x, y, z, u, v, w) + g:^H{x, y, z, u, v, w) = 0, 
d d d 

— A(a;, y, z, u, v, w) + ^^(^) ^) v, w) + ^^(^) ^) v.w) = 0. (22) 

This system of equation has a solutions depended from the arbitrary functions. 
As example the solution of the system depended from one arbitrary function lead to the six- 
dimensional space with the metrics 



W ^ I ^x{x)-^K{x) - (^^K{x)] I dx% 



i dvdw ^^^^ dudv ^ dudw ^ dv'^ ^ } ^^^^"^ 



+ {-^K{x)-f^K{x)-(^K{x)] \ dz'+ 
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\dudw dudv du^ dvdw (^)J ^^^^ 

-idxdu + dydv + dzdw (23) 

with arbitrary function K{x) — K{x, y, z, u, v, w). 
This metric is a Ricci-flat in a simplest case 

K{x, y, z, u, V, w) = F{x, y, z, u)w + v 

where the function F{x, y, z, u) is solution of the equation 

In more general case the solution of the system (22) lead to the space with the metric 

= (al^^ '''' + - ^''^^'^ - 

where x = {x,y, z,u,v,w) and L{x), M{x), N{x) are an arbitrary functions. 
The metric (24) has a six nonzero components of the Ricci tensor. 

The condition on the space (24) to be a Ricci-flat lead to the overdetermined system of equa- 
tions for the functions L{x), M{x), N(x) 

^xx -^xy ^xz 0? -^yy -^yz ^7 ^zz 

which is the six dimensional generalization of the four-dimensional Plebanski the second heavenly 
equation. 
In the case 

L{x, y, z, u, V, w) — a2{x, y, z)u'^v -\- a3{x, y, z)uv^ -\- a5 {x, y, z)u'^'w+ 

+ a()(x, z)uw + a8{x, y, z)w v + a9{x, y, z)wv , 
M{x, y, z, u, V, w) = b2{x, y, z)u^v + h3{x, y, z)uv^ + h5{x, y, z)u^w+ 

+b6{x, y, z)uw^ + b8{x, y, z)w'^v -\- b9{x, y, z)wv'^, 
N{x, y, z, u, V, w) — c2{x, y, z)u^v -\- c3{x, y, z)uv^ -\- c5 {x, y, z)u'^'w+ 

+c6{x, y, z)uw'^ -\- c8 (x, y, z)w'^v -\- c9{x, y, z)wv'^, 

where Cj, hi, Ci are an arbitrary functions we get the space which is the Riemann extension of 
corresponding affinely connected three-dimensional space ([5], [6]) 

ds^ = -2T^.dx'dx^'<bk + 2(ia;'^;. 
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One part of geodesic of such type of the space is defined by 



+ I (2 ag(,x,y, z) + 2 b5(,x,y, z)) ^y(s) + (-2 aS(x,y, z) + 2 c2(x,y, z)) -j-xis)] ^z(s)- 



(' 



ds ds J ds 

-2 a5 {x, y, z) j + (-2 b6 (x, y,z) + 2 c5 {x, y, z)) j ^vi^) = 0, 



^^,yis)-2b3ix,y,z)^-zis)j + 
+ (^{-2 c2{x, y,z) + 2 a3{x, y, z)) -^y{s) + (2 c3{x, y,z) + 2 b9{x, y, z)) ^x{s)j ^z{s)+ 
+ {-2b8{x,y,z) + 2a9{x,y,z)) (J^x{s)j ^y{s) - 2 c9{x,y, z) (J-^x{s)j =0, 

(^{-2 c5{x, y,z) + 2 b6(x, y, z)) ^^y{s) + (2 b8(x, y,z)-2 a9(x, y, z)) j ^^z(s)- 

-2a6(x,y,z) +{2 c6(x,y, z) + 2 a8(x,y, z)) ^y(s)-2 c8(x,y, z) (^j^x{s)^ + 

d^ , , 

Another part of geodesic is defined by the hncar system of equations 

+ A{x, y,z)-^ + B{x, y, z)^ = 

with some 3x3 matrix-functions A, B and ^ = {u,v,w). 

In considered case the conditions Rij — lead to the system of equations for the coefficients 

d 

Rxx = -4 a8{x,y,z)c3{x,y,z) - A —c8 {x,y, z) - 8 c8{x,y, z)a3{x,y, z)- 

—4 c6{x, y, z)b9{x, y,z)—A a8{x, y, z)b9{x, y,z) + 8 b8{x, y, z)a9{x, y,z) + 8 c9{x, y, z)c5{x, y, z)+ 

+8 c8{x, y, z)c2{x, y,z) — 4 c6{x, y, z)c3{x, y,z) — 8 c9{x, y, z)h6 (x, y, z) — 

-4 {a9{x, y, z)f - 4 {b8{x, y, z)f - 4 ^c9{x, y, z) = 0, 

Ryy = 8 a5{x, y, z)a9{x, y, z) — 8 a6{x, y, z)c2{x, y, z) — 8 a5{x, y, z)b8{x, y, z) — 
—4 a8{x, y, z)b5{x, y,z) — 4 a8{x, y, z)a2{x, y,z) + 8 a6{x, y, z)a3{x, y, z) — 
d d 2 2 

~^'g^"'^i^^y^^) ~ ^Q^^-^i^^y^^) ~ ^ i^^i^^y^^)) - ^ {c5{x,y,z)) - 
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—4 c6{x, y, z)a2{x, y,z) + 8 b6{x, y, z)c5{x, y,z) — A c6{x, y, z)b5{x, y, z) — 0, 
Rzz — 8 b2{x, y, z)b8{x, y,z) — 4 c3{x, y, z)b5{x, y,z) — 8 b2{x, y, z)a9{x, y, z) — 
—8 b3{x, y, z)c5{x, y,z) — 4 b9{x, y, z)b5{x, y,z) + 8 a3{x, y, z)c2{x, y, z) — 
—4 b9{x, y, z)a2{x, y, z) + 8 b3{x, y, z)b6{x, y,z) — 4 c3{x, y, z)a2{x, y, z) — 

-4 {a3{x, y, z)f - 4 ^h2{x, y, z) - 4 ^^^(x, y,z)-4 {c2{x, y, z)f = 0, 

Rxy — 4 c5(a;, y, z^a2{x^ y,z) + 4 b6{x, y, z)a9{x, y,z) — 4 b6{x, y, z)b8{x, y, z) + 

d d d 

+4 c5{x,y,z)b8{x,y,z) + 2—a9{x,y,z) + 2—a8{x,y,z) + 2 —c6{x,y, z)+ 

d 

+4 c8{x,y,z)b5{x,y,z) + 2—c5{x,y,z) + 4 a6{x,y, z)b9{x,y, z) - 8 c9{x,y, z)a5{x,y, z)- 

d d 
- c5{x,y,z)a9{x,y,z) - 2 —b6{x,y, z) + 4 a6{x,y, z)c3{x,y, z) -2—b8{x,y,z) = 0, 

d 

R^z ^ 4a3{x,y,z)b8{x,y,z) + 2—c2{x,y,z) + 4b3{x,y,z)a8{x,y,z)- 
—8 c8{x, y, z)b2{x, y,z)+4 c2{x, y, z)a9{x, y,z)+4 c9{x, y, z)a2{x, y,z)—4 c2{x, y, z)b8{x, y, z)+ 

d d 

+4 b3{x, y, z)c6{x, y,z)-4 a3{x, y, z)a9{x, y,z) + 2 ^^^(a;, y, z) + 2—b8{x, y, z)+ 

Odd 
"^^ a^^^'^^'^'^'' ~ 2— a,?(a;,|/,2;) -2—a9{x,y,z) + 4 c9{x,y, z)b5{x,y, z) = 0, 

Ryz = 4a5{x,y,z)b9{x,y,z) + 4b2{x,y, z)c6{x,y, z) + 4 a5{x,y, z)c3{x,y, z) — 
—4 a3{x, y, z)b6{x, y,z) + 4 b2{x, y, z)a8{x, y,z) — 4 c2{x, y, z)c5{x, y, z) — 

—8 a6{x, y, z)b3{x, y,z) + 4 a3{x, y, z)c5{x, y,z) — 2 —c5{x, y, z) + 

+2-^a3{x,y,z) + 2^a2{x,y,z) -2^c2{x,y,z) + 2^b5{x,y,z) + 
d 

+2—b6{x,y,z) + 4c2{x,y,z)b6{x,y,z) = 0, 
oz 

which can be of used for the extension of the Liouville theory of invariants of the second order 
ODE's on a three dimensional case. 

In particular the description of the projectively flat structures defined by the solutions of 
nonlinear p.d.e.'s connected with corresponding systems of ODE's is possible. 
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7 Eight-dimensional generalization 



By analogy the eight-dimensional generalization of the second Heavenly equation can be obtained. 
The metrics of corresponding space has the form 

/ g2 \ / g2 \ / Q2 Q2 92 \ 

+dxdp + rigriy + rfzrf-u + dtdv, (25) 

depending on the six arbitrary functions and x = {x, y, z, t,p, q, u, v). 

Remark that a Ricci tensor of the metrics (25) has a ten nonzero components. 

The generalization of the second Heavenly equations corresponds the system from a six equa- 
tions on the coefficients of the metrics (25). 
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